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We consider oscillators evolving subject to a periodic driving force that dynami¬ 
cally entangles them, and argue that this gives the linearized evolution around pe¬ 
riodic orbits in a general chaotic Hamiltonian dynamical system. We show that the 
entanglement entropy, after tracing over half of the oscillators, generically asymp¬ 
totes to linear growth at a rate given by the sum of the positive Lyapunov exponents 
of the system. These exponents give a classical entropy growth rate, in the sense of 
Kolmogorov, Sinai and Pesin. We also calculate the dependence of this entropy on 
linear mixtures of the oscillator Hilbert space factors, to investigate the dependence 
of the entanglement entropy on the choice of coarse-graining. We find that for almost 
all choices the asymptotic growth rate is the same. 
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I. INTRODUCTION 

The results of [1] state that, roughly speaking, if we have a quantum system with a 
finite-dimensional Hilbert space of states that we factorize into a product of Hilbert spaces, 

_ n = n A ®n B , (i) 
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then the typical pure state in H has very close to the maximal amount of entanglement 
allowed between Ha and Hb, and this is in turn maximized if dim(7-^ 4 ) = dim('Hs). We 
will call such factorization of H into “observable” and “non-observable” physics a coarse- 
graining of the system. This suggests that if we evolve a system randomly from an initial 
configuration with zero entanglement entropy, then it will eventually forget essentially all of 
the information of the initial state if we only measure observables sensitive to Ha- There are 
many studies of this kind of process in specific systems. The rate at which the entanglement 
grows towards saturation depends, in general, on the details of these systems, although some 
general bounds exist [2-5]. Linear growth in time appears in many systems, for example, 
studies of decoherence and quantum chaos [6-21], and quenches of extended systems [22-25]. 

When trying to apply these results in the context of black hole physics, we are usually 
confronted with two problems. First of all, the Hilbert space H is big. In the gauge/gravity 
duality [26] the dynamics takes place in an infinite-dimensional Hilbert space: it is the 
Hilbert space of a relativistic quantum held theory on the conformal boundary. 

A very naive application of the results of [1] would suggest that typical states have infinite 
entropy when splitting H in two pieces of the same size, since both are infinite dimensional. 
However, the notion of splitting along a random factorization has no meaning, because once 
we have factored into infinite dimensional pieces, we can factorize the pieces again: there is 
no natural notion of splitting in half. Thus, the question of the entanglement entropy for a 
typical state is ill-defined without additional structure on the Hilbert space. 

An example of such a structure is two operators algebras, one for Ha the other one for 
Hb- It is natural to do the splitting with respect to a choice of algebras with reasonable 
properties determined by features of the dynamics. Once that splitting is done, instead of 
computing the entanglement entropy of the typical state, we can compute the rate of growth 
of the entanglement entropy as governed by the dynamics and ask how this growth is affected 
by our choices of coarse-graining and dynamics. It is here that we need a model dynamics 
that is both tractable and generic. We have in mind two simple harmonic oscillators with 
two ladder operator algebras, and we will show in what sense a system like this can be 
considered generic. 

The second problem we find generically is that there is no obvious canonical splitting into 
two factors Ha and Hb, so one might expect that entanglement entropy based on some such 
splitting might depend substantially on the coarse-graining. If we also define the scrambling 
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rate as the slope of the entropy growth, one might worry that there is no objective way 
to quantify it. This would make it very hard to understand in what sense black holes are 
fast scramblers [27, 28], when we think of the evolution in terms of a dual quantum held 
theory. In this sense, it is natural to ask if there is a universal result where the details of 
the factorization do not matter too much. Our main motivation is to eventually formulate 
the fast scrambling conjecture on a rigorous footing, but to do so, we need to be able to 
apply the methods that could characterize scrambling to fairly generic dynamical systems 
to which we associate an infinite-dimensional Hilbert space. 

The purpose of this paper is to analyze the scrambling rate, i.e., the entropy growth rate, 
in a toy model of chaotic dynamics that iterates a relatively simple unitary evolution operator 
on an infinite-dimensional Hilbert space %. The Hilbert space will be further decomposed 
into a product of two infinite-dimensional Hilbert spaces. This is done for a closed system, 
and we will study the dependence on the choice of coarse-graining and initial state. The 
idea is to study the quantum counterparts of closed Hamiltonian chaotic dynamical systems 
with finitely many degrees of freedom, similarly to previous studies of closed systems in the 
quantum chaos and decoherence literature [9, 12, 18]. 

So long as these classical systems have bounded trajectories (for example if the regions 
with bounded energy have finite volume), then they are expected to have a dense set of 
periodic trajectories. We assume that these periodic trajectories encode all the important 
information of the dynamical system, and that any classical initial condition is sufficiently 
close to such a periodic orbit, in line with periodic-orbit theory in quantum chaos [29-31]. 
The evolution of the classical system for such an initial condition can then be understood, for 
some time, by the linearized perturbations around the corresponding periodic orbit. We can 
ask how these perturbations grow in time and estimate the Lyapunov exponents of the full 
system from such an analysis. Given such Lyapunov exponents, one then has the classical 
entropy growth rate, essentially the Kolmogorov-Sinai (KS) entropy. 1 

The main point of this paper is that the entropy growth rate is the same in the classical 
and quantum systems and that it is essentially independent of the coarse-graining and initial 
state (assuming that the latter is a sufficiently good approximation to a classical state). Our 
results are consistent with the seminal work on open systems by Zurek and Paz [6], where a 
1 See [31-33] for general discussions on chaotic Hamiltonian dynamics and more precise definitions. The 
technical definition of the KS entropy requires a finite measure space, while we are working with a phase 
space and Liouville measure that are only c-finite, but our results do not depend on this issue. 
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preferred factorization into a degree of freedom and an environment is given, but note that 
we explicitly find that the entanglement entropy growth rate grows as the sum of the positive 
Lyapunov exponents, not just as the largest one. The original work of Zurek and Paz [6] 
had only one such positive exponent, and similarly with [12, 18], though other studies have 
also found growth rates equal to the full sum [7, 9] . We are also consistent with other earlier 
studies of entanglement entropy growth, but distinguished from them in that our system is 
not coupled to an external environment [6-8, 10, 14, 19, 21], our results are not perturbative 
[13] or numerical [34, 35] and we do not work with a finite-dimensional ffilbert space [9, 11] 
or invoke a random matrix or semi-classical approximation [15-17, 20, 36]. 

We also go further in that we explicitly show that the choice of coarse graining in this 
closed system almost always matters very little. We consider all factorizations of the ffilbert 
space in which each factor corresponds to half the degrees of freedom. 2 This is the closest 
analog we can imagine in this class of systems to saying that Ha and He are the same size 
(since both are infinite dimensional), while ensuring that the observables in Ha commute 
with the observables in He- We focus on the minimal system for which this splitting can 
be done, namely, two degrees of freedom, but the methods used here can be generalized to 
more degrees of freedom. 

II. A SIMPLE STROBOSCOPIC DYNAMICS 

As described in the introduction, we are interested in the entanglement entropy growth 
for a quantum system with an infinite-dimensional Hilbert space after acting repeatedly on 
it with a unitary operator on a preferred initial pure state. The quantum system is supposed 
to arise from quantization of a classical (chaotic) dynamical system. We consider an initial 
pure state as close to a classical initial condition as possible given the uncertainty principle, 
like a coherent state of a harmonic oscillator, since this should allow us to clearly see the 
relationship between classical and quantum entropy growth. 

We are not able to solve this problem in general, so we solve a highly simplified model as a 
first step. In chaotic Hamiltonian dynamical systems with bounded phase space trajectories, 
e.g., those whose constant-energy surfaces are compact, there are many periodic solutions 
2 More precisely, over a factor spanned by the states generated by the quantization of a single particle 
coordinate q and its canonical conjugate p, which we combine into ladder operators q ± ip. 
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(in highly chaotic situations they are dense in the phase space), and we will choose our 
initial condition to lay very close to one such periodic trajectory, with period T, but in a 
chaotic domain of phase space, i.e., away from invariant (Kolmogorov-Arnol’d-Moser) tori. 
If we choose such an initial condition and wait a time T, we will end up close to our initial 
condition and we have a canonical transformation that has a fixed point set: the periodic 
trajectory itself. Other nearby periodic trajectories will have much longer periods in general, 
so we can ignore them if we only study the system at times comparable to T or perhaps one 
order of magnitude larger. This is all we need for most of our results to hold, but see the 
end of Sec. IV for a more quantitative discussion. 

The simplification of the dynamics occurs by considering an infinitesimal neighborhood 
about some point on this periodic orbit, and assuming that this can be treated with a lin¬ 
earized, tangent-space approximation to the dynamics. This approximation should keep all 
the essential features of the dynamics in an expanded, finite neighborhood, by the Hartman- 
Grobman linearization theorem. 

Note that by translating the solution by a small amount of time again leads to an exactly 
periodic solution with period T, so there is no stretching in this direction. Geometrically, 
this is the direction along the periodic orbit that pierces the small neighborhood we are 
considering. This direction that experiences no growth leads to one vanishing Lyapunov 
exponent. There is a second direction that experiences no growth, which is related to 
conservation of energy. These two form a conjugate pair. We will consider only directions 
linearly independent of these directions, and so will generally have an even number of non¬ 
zero Lyapunov exponents. 

Thus, we are led to consider a simplified problem where we have a linear canonical 
transformation in phase space such that the origin is fixed. Eventually we want to treat the 
system quantum mechanically, and we will form a quantization of this linear phase space 
in terms of operators q s ,p s that transform linearly into each other. We do this by first 
choosing a set of Darboux coordinates, q s ,p s , centered on the fixed point, that parametrize 
the tangent space of the fixed point and simplify the symplectic form to A dp s . 

Then we can associate a copy of the ladder operator algebra with each pair of canonical 
coordinates: 

a s = (q s + ip s )/V2 and a\ = (q s - ip s )/V2 . (2) 

Evolution in our model is effected by a canonical linear transformation, which can be 
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represented by a real matrix that preserves hi. In matrix notation, preserving hi corresponds 
to the equation 

M t QM = Q, (3) 

where we choose to represent hi as an antisymmetric matrix and M T is the transpose of M. 
We can then decompose the linearized phase space into subspaces that diagonalize M. 

Assuming that we have an eigenvalue A, since M is real, we must have that A is also an 
eigenvalue. Furthermore, preserving can also be written as 

M t = f, (4) 

which shows that M T is conjugate to M~ l . Since the transpose of a matrix and the matrix 
have the same set of eigenvalues, any eigenvalue of M T corresponds to an eigenvalue of 
M. The eigenvalues of M T are the inverses of the eigenvalues of M. We thus find that to 
any eigenvalue, we can associate a four-dimensional subspace of phase space, the (possibly 
complexified) linear space spanned by the eigenvectors with eigenvalues A, A, A -1 , W 1 , which 
are generically distinct. The eigenvalues cannot be zero because is non-degenerate. 

The action of is closed on these linear subspaces. Consider two vectors vi,v 2 that are 
eigenvectors of M. Then we have that 

= Aiujfhq = V2(M~ 1 ) T £lv\ = (M~ l v 2 ) T Vtvi = X^v^Qvi , (5) 

so if the matrix element of ujhhq is non-zero, we find that A 2 = A]" 1 . That is, we find 
that the linearized phase space can generically be decomposed into four-dimensional linear 
subspaces, such that the action of M and closes on these subspaces, so long as there are 
no degeneracies in the eigenvalues. 

There are three different cases of degeneracy to consider. First, assume that A = A. 
Then the eigenvalues are real, and the subspace is only two dimensional. These will be 
called unstable directions, because upon iteration the distance from the origin will grow 
indefinitely. The second possibility is that A = A -1 , which also gives a two dimensional slice 
of phase space. Then the eigenvalues lie on the unit circle, and the direction is said to be 
stable. Finally, assume Ai = Aj" 1 , in which case A = ±1. These necessarily come in pairs of 
eigenvalues because is antisymmetric and non-degenerate, so there must exist a vector v 2 
to any eigenvector v\ such that Affhq 7 ^ 0. By the argument above, if v 2 is an eigenvector 
of M itself, it has the same eigenvalue. 
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This third case lies at the intersection of stable and unstable subspaces. These will be 
called marginally unstable. It is in this case that more exotic things can happen because 
the matrix M might not be diagonalizable, though it can always be put in upper-triangular 
form. A simple example is a two-dimensional phase space where we get the transformation 



which has two identical eigenvalues, but the matrix does not have two eigenvectors. Itera¬ 
tions of the transformation lead to a constant drift 

q n = qo + anp (7) 

depending on the conjugate momentum. This corresponds to sub-exponential growth of the 
distance between initially nearby trajectories. These blocks are of vanishing measure in the 
set of all possible linear canonical transformations, but they might appear as non-trivial 
limits. It is an example like this that corresponds to the pair of coordinates consisting of 
energy plus the direction of flow along the periodic orbit. The role of p is played by the 
energy, while go is the linearized coordinate along the flow. 

This is also relevant for integrable systems with bounded orbits, which are characterized 
by having all directions stable or marginally unstable. For the rest of this paper we work in 
a generic case with no degeneracy. 

III. QUANTUM LINEAR DYNAMICS FOR TWO DEGREES OF FREEDOM 

WITH A CONSERVED U(l) CHARGE. 

The entanglement entropy is calculated by first factorizing the infinite-dimensional 
Hilbert space, generated by the above ladder operators (see Eq. (2)), into two infinite¬ 
dimensional factors. We are effectively treating the q s and p s operators as position and 
momentum operators of harmonic oscillators. We will trace over one factor to compute the 
entanglement entropy. The idea is to show that the classical entropy rate (the sum of the 
positive Lyapunov exponents) is equal to the growth rate of the entanglement entropy for 
late times, for almost all choices of coarse-graining (factorization). 

To simplify matters further, we will require that the phase space be only of dimension 
four. That is, we consider the phase space to be described by two q,p pairs of canonical 


coordinates, or in the quantum theory, by two harmonic oscillator algebras. As a further 
simplification, we will impose a U(l) symmetry, so that the linear action preserves a U(l) 
charge, with the two raising operators having opposite charge ." 5 This choice is partly from 
our desire to understand entropy production near black hole horizons a la Hawking [37], 
where there is a distinction between particles and antiparticles in terms of modes that fall 
into or escape from a black hole. The ones that fall into the black hole are considered 
unobservable and one needs to trace over them. In the end, the main advantage of the U(l) 
symmetry is that it leads to simple but illustrative computations. 

Furthermore, this requirement of the U(l) symmetry is not as constraining as it might 
seem. If we assign a charge 1 to the coordinate whose eigenvalues under the evolution matrix 
M is A, we can also assign charge +1 to the coordinate associated to the eigenvalue A -1 , 
and charges —1 to the other two coordinates, associated with the eigenvalues A and A -1 . 
This charge is preserved by the canonical transformation. Thus, a conserved charge can be 
defined whenever A is a generic complex number. 

Let us now write the most general form of such a linear dynamics in the quantum theory, 
using two ladder operator algebras, a, gA, 6 , b\ with both gA and b of charge 1 , and working in 
the Heisenberg picture of quantum dynamics. The most general canonical transformation 
will take gA to a linear combination of gA and b, and similarly, b will be taken to a linear 
combination of b and gA . For example, 

gA = Aa^ New + Bb New . (8) 


The most general such transformation can be parametrized as follows (see Appendix A for 
the derivation) 




= exp (i/3) 


e ld cosh p e 1 ^ sinli p 
e~^ sinh p e~ l6 cosh p 



( 9 ) 


The transformation of a, fA follows from taking hermitian conjugates. This is a Bogolubov 
transformation. 

This shows that the associated group is a U(l) x SL(2, M), where the U(l) corresponds to 
the phase exp (i/3). Indeed, if we forget about /3, the matrix in Eq. (9) has determinant equal 


3 This will effectively reduce the possible evolution operators to consider from all of Sp(4,M) (this includes 
consistency with the adjoint operation) to the subgroup SL(2,R), as we see in more detail below. 
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to one, and trace given by t = 2 cosh p cos 9. The eigenvalues of that matrix are solutions of 

A 2 — f A + 1 = 0 . (10) 

There are two cases to consider: \t\ > 2 and \t\ <2. In the hrst case the eigenvalues are real 
and inverses of each other, so we parametrize them as ±exp(p) and ±exp(—p), whereas in 
the other case the eigenvalues are on the unit circle and are exp(±?’0) . Therefore we write 
t = ±2 cosh p or t = 2 cos 6 for both cases. 

For the full matrix we need to multiply both eigenvalues by exp (i/3), thus we can make 
sure that t > 0 by changing the sign of exp (i/3) if necessary. When the eigenvalues are on 
the unit circle, the system is said to be stable. We are mostly interested in understanding 
the unstable case, but we will work in the general setup anyhow. 

The next thing we need to do is to choose an initial state. This will be the vacuum of 
the a, b oscillators, characterized uniquely by 

a \ 0) 0 = 6 | 0) 0 = 0 . ( 11 ) 

We want to know what the new state will be after application of the canonical transforma¬ 
tion. The answer follows from expressing 

6|0)o = (-D&New + C a New)|0)o = 0 • (12) 

The state can be interpreted as an “eigenvector” of the lowering operator & New , with an 
operator valued eigenvalue proportional to ajj ew , so it can be written in a similar form as we 
do with coherent states, namely 

|0 ) 0 oc exp(—(7.D - 1 al&l)|0)i (13) 

and in general 

|0 ) 0 oc exp(—C'„D“ 1 a t 6 t )|0) n (14) 

Notice that this expression is independent of /3, so /3 plays no significant role in the rest of 
the paper. 

Our notation is as follows. The vacua |0)„, will denote the vacuum of the oscillator algebra 
a, a', 6 , /A after n iterations of the unitary evolution characterized by the matrix M. The 
initial state | 0 )o is a squeezed state in terms of the new oscillator algebra after n iterations. 
The preserved charge is 


Q = a'a — b^b 


( 15 ) 
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we will say that a t creates particles and b1 creates antiparticles just as is done in conventional 
quantum field theory. 

To choose different initial states, we do another Bogolubov transformation to a different 
set of raising and lowering operators at time zero, such that the Bogolubov transformation 
preserves the charge. This will change the matrix M by conjugation, but the eigenvalues of 
M will stay the same. The new state will be the vacuum of the new set of oscillators, which 
can be interpreted as an initial state which is a squeezed state in terms of the old algebra. 
All of the states we have chosen as initial states are minimum uncertainty states centered 
at the fixed point. 

Notice that if we fix t = 2 cosh p, there is a particularly simple matrix with the right 
characteristic polynomial. This matrix is 


M = 


(16) 


A B\ ( cosh p sinli p 
C DJ \ sinh p cosh p* 

Any Bogolubov transformation with the same t is actually conjugate to this matrix. The 
powers of M are given by 


M n = 


A n B n 
C n D n 


A B\ ( cosh np sinh up 
CD) \ sinli np cosh np 


(17) 


so that 

C n D~ l = tanh np . (18) 

Similarly, for \t\ < 2, t — 2 cos 9, a simple M with the right characteristic polynomial is 
given by 

exp {id) 0 


M = 


A B 
C D 


0 exp (—id) 


(19) 


so in this case the oscillator algebras don’t really mix. 


IV. ENTROPY GROWTH I: TRACING OVER ANTI-PARTICLES 

Our goal is to analyze the entanglement entropy of a squeezed state as described in 
equation Eq. (13). This is, we want to understand the entanglement entropy of a state of 
the form 


\Q) oc exp(cra^)|0) 


(20) 
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where a is a complex number. To do that, we need to choose a factorization into a product 
of two Hilbert spaces. Since in this setup it is natural to separate the Hilbert space according 
to the oscillator algebras of a, b separately, here we will trace over the b oscillators. This 
is what one would usually do in the case of black hole physics. Then we will compute the 
growth of the entropy as a function of time according to the iteration in Eq. (14). 

The state | O) needs to be normalized. We will consider normalized states | n,m) for the 
two sets of oscillators. In terms of these, the states a' k W k |0) = k\\k, k) have norm k\. With 
these we find that 

k 

exp(cra t & t )| 0 ) = = ^^k.k) . ( 21 ) 

The norm of the state is therefore 

I exp(aaW)|0)| 2 = V \a\ 2k = - \ (22) 

1 — \ar 

k 

and it is finite only if \a\ < 1 . 

When we trace over the b degrees of freedom we get a reduced density matrix for the a 
degrees of freedom given by 

p a = (l~\a\ 2 )^2\a\ 2k \k)(k\ . (23) 

This is the same as a thermal density matrix for a harmonic oscillator a, a' with the Boltz¬ 
mann factor identified as exp(— /3hw) = |a| 2 . So from the perspective of the a oscillator, the 
dynamics is pumping in heat. 

The entropy of this density matrix is the entanglement entropy of the subsystem of the 
a oscillators and it is then given by 

S = — log(l — x) -—— logrc (24) 

x — 1 

with x = | ck| 2 . Notice that this entropy is independent of the phase of a. 

We now apply this result to the squeezed state that arises from the iteration of our unitary 
dynamics in the special case of Eq. (17), which results in the identification 

\a n \ = tanhn/3 (25) 


so that for large n we have 

x = tanli 2 np 


1 - 


cosh(np)" 


->• 1 - 


exp( 2 np) 


( 26 ) 
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and we find that 

S n —» log(exp(2np)/4) + 0(1) = 2np + 0(1) . 

We thus find that the asymptotic entropy growth per unit time is characterized by 


AS 
A t 


1 AS 


2 p 
T 


(27) 


(28) 


T An 

where we have inserted the time scale of the periodic trajectory alluded to in the previous 
section. 

This shows that the entropy growth is asymptotically linear as a function of time, and 
that the growth is controlled by the eigenvalues of the matrix M. We will later interpret 
this result in terms of Lyapunov exponents. Our goal is now to show that this result is 
essentially independent of the choice of basis for the harmonic oscillators. The first step 
is to consider a more general matrix as in Eq. (9). Since t = 2 cosh p = 2 cosh p cos 9, we 
immediately find that p > p and that therefore the entropy for other choices of initial states, 
related to the vacuum by the action of such a general matrix, is larger than that for the 
special case chosen above. We will find that this only affects the 0(1) piece, not the leading 
n-dependent piece. 

Since each such more general matrix is conjugate to M, we can write its nth power as 
follows: 

e *6»n gQgp p n e *<+ gink p r 


sinh p n 




cosh p n 


= TM n T - 1 


with 


T = 


e* 72 cosh 7! e* 73 sinh 

e - * 73 sinh 71 e - * 72 cosh 71 

It follows from a straightforward computation that 


(29) 


(30) 


coshp n exp(f6 ) r i) = cosh(np) + cosh^) sinh(7 1 )(exp(i7 3 — i^ 2 ) — exp(ry 2 — *73)) sinh(np) 
= cosh(np) + 2 i cosh^) sinh^) sin(7 3 — 72) sinh(np) . (31) 

With this information we find that | coshp n | 2 is modified asymptotically as 

12 1 + | sinh( 2 yi) sin (73 - 72 ) | L> 


cosh p n | ~ 


exp(2np) 


(32) 


and one can show that in the formula determining the entropy this only affects the order 
0(1) terms by log(l + | smli^yx) sin(7 3 — 72)| 2 )- This is independent of n, so we still have 
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the result 

S n = 2np + 0(1), (33) 

independent of the choice of initial state, since such a choice amounts to another conjugation, 
as discussed below Eq. (15). Out of these choices of initial state and evolution operator, 
the entropy is smallest for the special matrix M we chose, but asymptotically this difference 
in negligible. We see here that the n at which the linear growth term dominates depends 
on the size of the y, and p. In practice, if the y* and p are all 0(1), then even for n < 10 
the linear approximation is quite good, justifying our remarks in Sec. II that we can ignore 
other potentially nearby periodic orbits. If one considers more extreme values, the value of 
n at which we reach the asymptotic behavior above depends on details. 


V. ENTROPY GROWTH II: THE GENERAL CASE 


One can imagine that, in general, different factorizations of the Hilbert space will lead to 
different answers for the entanglement entropy and that this might affect the entropy growth 
we computed in the previous section. The purpose of this section is to show that this is 
not the case: if one uses a different linear combination of raising operator states to define 
particles versus antiparticles, we will see that in general they all have the same asymptotic 
growth of the entropy. We will also generalize the results to include general canonical linear 
transformations, not only those that preserve a U(l) charge. 

So far we have chosen to trace over the b oscillators for a pure state of the form 

10) cx exp(craW)|0) (34) 


to compute the entanglement entropy growth. The simplicity of this particular choice is 
that the associated density matrix for the a oscillators is already diagonal in the number 
basis when we trace over the b oscillators. We will now study a more general case where 
this diagonal form is not readily available. 

The idea here is to rotate the oscillator algebra as follows 



cos 9 sin 9 
— sin 9 cos 9 



(35) 


Here the subscript r is reminding us that the basis is rotated. Such a rotation is also 
a Bogolubov transformation that preserves the ladder operator algebra, but it does not 
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correspond to an evolution like the one we have considered before, since it violates the U(l) 
symmetry, i.e., it does not preserve the Noether charge Q = a)a — tfb (except for 9 = 0). 
The vacuum |0) does not change under such a rotation, but the new way of writing the state 
is more complicated. It can be thought of as a redefinition of particles and anti-particles at 
a given time, or a different coarse graining of the system into object and environment. The 
underlying state can be written 

|fi) = C exp [a(— cos 9 sin da\ 2 + (cos 2 6 — sin 2 9)albl + cos 6 sin 6bl 2 )] |0) . (36) 

We now want to trace over the b r oscillators and compute the density matrix for the a r 
oscillators. We can also consider a more general SU(2) rotation for equation (35). This will 
modify various relative phases in equation (36), but the structural form will be very similar, 
since the additional phases can be eliminated by changing the phase of the and the 
aj., bl operators. So, we can assume equation (35) without loss of generality. 

A convenient way to compute the reduced density operator is to use the holomorphic 
representation of the oscillator Hilbert space (also called the Bargmann representation). 
Here ~ z and a ~ d z . We can arrive at this representation by noticing that for any 
entire function i^{z) there is a corresponding state \ij}) = ^(a')\0) and that the annihilation 
operator acts as d/dz, i.e., a\^) = -0'(at)|O). For two states |</>) and |-0) we can write their 
inner product as 

(4>\^) = — J 4>(z)ip(z)e~ zz dzdz. (37) 

See Appendix B for further details. Thus we may represent our state by 

10) —> fi(£, w) = C exp [a(— cos 9 sin 9z 2 + (cos 2 9 — sin 2 9)zw + cos 9 sin 9w 2 )] (38) 

and similarly 

(Q —>■ f \{z,w) = C' ex p [a(— cos 6* sin #7 2 + (cos 2 9 — sin 2 9)zw + cos 9 sin 9w 2 )\. (39) 

Then we can implement the trace over the b oscillators via the integral 

p a = Tr b |Q) (0| —» p a (z,z) — — J Q(z,w)Cl(z,w)e~ ww dwdw, (40) 

with an abuse of notation by writing the holomorphic representation of p a as p a (z, z). Since 
the integrand is the exponential of a quadratic polynomial, we can perform the integral. For 
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simplicity of notation let us assume that a € M, which can be achieved by changing the 
phases of the variables z, w. 

We get 


Pa(z,z) = 


1 — or 


\Jl — oi 2 sin 2 29 


exp 


a- 


(a 2 - 1) sin(20)(^ 2 + z 2 ) + 2acos 2 (29)zz 


(41) 


2(1 - a 2 sin 2 29) 

where we used \C'\ 2 = 1 — a 2 (see Eq. (22)). To get the ( m,n ) matrix element of p a we 
simply read off the coefficient of z m z n in the above, times \/m\n\. It is useful to rewrite p a 
as 

p a (z,z) = B oe B ^ z2+ ^ +B2Z * , (42) 


where 


Bo — , , 

V 1 — a 2 sin 2 29 

a(a 2 — 1) sin 20 

1 2(1 — a 2 sin 2 29) 

a 2 cos 2 29 

2 1 — a 2 sin 2 29 


(43) 

(44) 

(45) 


This is the same form of density matrix as that obtained in [38], where they were able to 
compute the entanglement entropy. We reproduce their computation in detail in Appendix 
B and fold 


S = — Tr (p a log p a ) = (x + I) log(x + I) - (x - |) l°g(x - \) , (46) 

where _ 

a 2 cos 2 29 ^ 1 
(1 - a 2 ) 2 + 4 ' 

What is important for us is that for large times we have that a n —> 1, so substituting a n for 
a in Eq. (47), the resulting value of \ n becomes large and is dominated by the singularity in 
the denominator inside the square root. For large y, we can expand the entropy as follows 




S ~ (x + \) lo g(x) “ {X |) log(y) + 0(1) « log(x) + 0(1) (48) 

and log(x) ~ log(l/(a 2 —1))+0(1). Substituting \a n \ = tanh np we get the same asymptotic 
answer as before, but the terms of order one now act to lower the entropy rather than to 
raise it. See Fig. 1 for a plot of the entropy as a function of n for various choices of 9. 
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Figure 1. The von Neumann entropy of p a , as given by Eq. (46), with a n = tanhn, as a function 
of n, for 8 = 0 , (0.5)7r/4, (0.9)7 t/4 , (0.99)7 t/4 (topmost to bottom-most curve) . 


It should be noted that this reasoning is correct so long as \A 2 \ ^ 0, that is, cos2 6 ^ 0. 
For that very special case, the coefficient of aJtf in the squeezed state in Eq. (36) vanishes, 
and the state is a product state between the two Hilbert spaces. In that case the entropy 
vanishes and x — Thus, there is one special selection of coarse graining that does not 
produce entanglement entropy. 

The other case we have not treated so far is when instead of having a complex eigenvalue 
for the canonical transformation (and its three images), we have two different real eigenvalues 
and their inverses. In such a case, to each pair of eigenvalue we can associate a q,p pair in 
phase space, and ladder operators cj, c i: as before. The Bogolubov transformations will mix 
the c, c f operators for each i independent from each other. Such a transformation can be 
parametrized as 

c = (cosh p)c New - (sinh p)cjj ew (49) 


and this leads us to a squeezed state given by 


|0) o oc exp ( f (tanhpjJcJ 2 + Utanh p 2 )<% 


(50) 


This is a product state, and there is apparently no entropy production. In the holomorphic 
basis, this corresponds to a state 


|0) 0 oc exp ( -(tanhp^^ 2 + -(tanhp 2 )^ 2 


(51) 
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and to eigenvalues exp(±pi i 2 ) for the different directions in phase space. 

Upon iteration, we would get 

= exp (|x 2 + |V) . (52) 

What is important to notice is that so long as p 12 are non-zero, the large iteration limit 
produces coefficients that tend to 1/2, and we have introduced si,s 2 in order to write 
compact expressions later. 

Just like before, we can rotate the algebra and ask what the dependence of the entropy 
growth on the choice of rotation parameter is. The rotations that are available to us given 
c\ 2 belong to SU(2). We are interested in the large iteration behavior. However, notice that 
if we do a linear transformation between c\,c\ with real parameters, we get a cancellation 
of the cross term in the squeezed state, producing a state that is still factorized. Thus, it is 
more convenient to consider a complex transformation for illustration purposes: 


cos(0)^i + i sin(0)wi 

(53) 

isin(0)^! + cos(0)tci 

(54) 


Alternatively, one can take p 2 — P 2 and use a real rotation between z,w to go to the 

new basis. Integrating out the new variable w 1 again produces a density matrix of Gaussian 
form. The corresponding squeezed state in the new basis is 

|0)o oc exp (cos(0)zi + i sin(0)tci) 2 + -^-(i sin(0)zi + cos(6 l )tci) 2 j (55) 

Call the quadratic form in the exponent Q(z\,Wi). To obtain the density matrix for z±, we 
take the partial trace over the w\ modes. This is given by 

Pa{zii Al) — y'd 2 Wiexp[Q(^i,«; 1 ) + Q{z 1 ,w 1 ) - w^] . (56) 

The quantity in the exponential is a gaussian with a shift in wi,Wi. When we integrate over 
wi,wi, we replace them by the values at the critical point of the Gaussian for fixed Z\,Z\. 
The final answer for this quadratic form is not particularly illuminating. Instead, we quote 
directly the value of y 2 which is 

2 (1 + S 1 S 2) 2 — (®i + ^ 2 )^ cos 2 (20) 

r= 4 ( 1 -*;)(!-»!) 


|0) 0 oc exp ( -(tanhnpi)^ 2 H—(tanh np 2 )v/ 


(57) 
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Notice that this is symmetric in the exchange ,s'i -h- s 2 , as it should be: the entanglement 
entropy of subsystem a is the same as the entanglement entropy of subsystem b. Inserting 
Si = tanh(npi) and s 2 = tarih(n/> 2 ) we get that 

X 2 = yx[3 + cosh(2n(pi + p 2 )) - 2cos(46») sinh(ra(pi + p 2 ) 2 )] , (58) 

lo 

and as a consistency check we find that x 2 |e=o = 1/4 (for which no entropy is generated). 
Taking the large iteration limit, we see that the entropy S ~ log(x) generically grows as 

S = logexp(n(pi + p 2 )) + 0(1) ~ n{p i + p 2 ) (59) 

except on a set of measure zero (where 008(46*) = 1). Notice also that when pi = p 2 = p 
we recover the results of the previous section. In general, the entropy grows as the sum 
of the logarithm of the two large eigenvalues of the linear symplectic transformation. This 
result also includes the case where one pair of eigenvalues is unitary and the other is real. 
The unitary (stable) eigenvalue pair does not contribute to the entropy because in that case 
s 2 = 0. 

As a further comment, we should note that all of the calculations work essentially the 
same way if we start initially with a coherent state, rather than just a vacuum. This is 
because a coherent state equation a\f) = f\f) can be thought of as defining a new set of 
ladder operators via a shift a = a — f and a) = a' — /*. This is an automorphism of 
the algebra of raising and lowering operators (it preserves the commutation relations and 
the adjoint operation), so it can be used to show that any such coherent state is like a 
vacuum, and any “shifted” squeezed state is a squeezed state with respect to a new set of 
shifted operators. Essentially, if we extend the linear problem of Bogolubov rotations to 
allow these translations, we get an affine group. As such, any element can be composed of 
a “rotation” (a symplectic matrix) plus translation and the translation part can be moved 
from multiplying on the right to the left without affecting the rotation. This is, if we think 
of / as a translate of the vacuum 

l/}-T(/)|0> (60) 

by a unitary operator T(/), and then apply a Bogolubov rotation R{0) to it (as an active 
transformation), we find that since 


R(0)T(f) ~ T(f)R(0) 


(61) 
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we get that 


R(9)T(f)\0)~T(f)R(9)\0) 


(62) 


and the effects of T(f') can be undone on the algebra of operators of the target by a passive 
transformation (it acts on the definition of the algebra, not the states). These translations 
furthermore preserve the splitting of the Hilbert space that we have described so far (for 
any choice), so they do not entangle the Hilbert spaces further and when we compute the 
entanglement entropy, it doesn’t affect the answer. 


VI. COMPARISON WITH THE CLASSICAL DYNAMICS 


Consider a general chaotic Hamiltonian dynamical system with finitely many degrees of 
freedom, corresponding to canonical pairs qi,Pi , and with a time independent Hamiltonian. 
Assume that we have a periodic trajectory with period T, characterized by q[°\t ), p\°\t) and 
consider the evolution of infinitesimal deviations, 5qi(t),5pi(t), from that trajectory. The 
equations of motion are given by 


Qi 

Pi 


d H 

dpi 
d H 
dqi 


Using qi = q + 5qi, p. t = pf^ + dpi, we can Taylor expand 


(63) 

(64) 


• • (o) I r • dH 

qi = q\ +dqi = — 

dpi 


£ 


dH 


dpidsj 


q(0) ,p(0) 


5sj + 0(5s 2 


(65) 


lg(°),p(0) 3 

and similarly for p. Here the sum Sj is over all the canonical variables. Since q ^ is a 
solution of the equations of motion, we find that the evolution of the 5sj is given by a first 
order set of coupled homogeneous linear differential equations. These have time dependent 
coefficients characterized by the quantities 




± 


OH 

dsidsj 


q(0),p( o) 


( 66 ) 


which are also periodic with period T. There is a monodromy matrix $(T) associated 
with this problem, a Floquet problem, that expresses any solution after time T in terms of 
known solutions with t < T, in the form u(t + T) = <3>(T) • u(t). Here u(t) is a vector of 
linearly independent solutions to the equations of motion, and <f>(T) is a representation of the 
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Floquet operator (see, e.g., [39] for a review of this formalism). This can be iterated, so that 
u(t + nT) = $(T) n • u(t). The monodromy matrix $(T) plays the role of M (the evolution 
matrix discussed in section II). One of the eigenvalues of <E>(T) is 1, corresponding to the 
the eigenvector representing the infinitesimal deviation Ssi = s\°^8t. This just corresponds 
to changing t —> t + 5t. 

As above, we assume that we are in a generic case, so that we can decompose our vector 
space along eigenvectors of $(T) and only keep the span of the vectors that correspond to 
eigenvalues different than 1. This is the space which we assumed is four-dimensional in the 
previous sections. 

One can do something similar on Hamiltonian systems with a periodic driving and with 
a periodic trajectory, which are associated with a time-dependent Hamiltonian H(t). Such 
systems can be extended from N to N + 1 degrees of freedom, where one takes H^ ew = 
p s + H(s) with (s,p s ) the new pair of canonical variables. The equations of motion of s 
are s — 1, so s is essentially time and now any such time dependence is due to solving the 
evolution of a time independent system. These systems are said to have N + 1/2 degrees 
of freedom [32, 33]. They have the advantage that eigenvalue that is equal to unity can 
be associated with the pair ( s,p s ) itself, so the associated monodromy matrix $(T) can be 
thought of as non-degenerate and even-dimensional from the get go. 

We now consider the Lyapunov exponents for this system. Let us assume we are in the 
generic case, and write the eigenvalues of <f>(T) in groups of four, A, A, A -1 , A -1 , all distinct, as 
above, with A = exp (p+i/3). We take p > 0 since that is the case in which we are most inter¬ 
ested. Let us further work in a (possibly complex) basis that diagonalizes <f>(T), so that the 
Jacobian matrix of the dynamics in this basis is simply J = diag(e p+l/3 , e p ~ l/3 , e~ p ~ l/3 , e ~ p+l /3 ). 4 
The Lyapunov exponents are the logarithms of the eigenvalues of the matrix 

L — lim ((J n ) t j 71 ) 1 / 2 " — diag(e p , e p , e~ p , e~ p ) , (67) 

n—>oo v 7 

so for each such group one has the Lyapunov exponents (p,p, — p, —p). The exponents for 
the continuous-time system have to be scaled by T, so this indicates the rate of separation 
of nearby points grows exponentially in time as exp (pt/T). 

For sufficiently uniform chaotic systems on domains with finite measure the Kolmogorov- 
Sinai (KS) entropy is equal to the sum of the positive Lyapunov exponents K t (times their 


4 We are effectively using covariant Lyapunov vectors [40, 41] (see also [31] Ch. 6). 
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degeneracy d,): 

/Jks = ^2 diKi ’ (68) 

K ;>0 

by Pesin’s theorem. 5 This is the sum of all expanding Lyapunov exponents. Roughly 
speaking, the KS entropy expresses the rate at which we lose information about the initial 
conditions due to the chaotic dynamics. One starts by subdividing the set of possible initial 
conditions into a coarse-grained set of small volume elements, or cells. The KS entropy 
computes how many such cells are covered by our initial volume when we evolve it to 
asymptotically large times and take the size scale of the coarse-graining to zero, see, e.g., 
[32, 44], 

We can use the same intuition in our case, even though we are considering dynamics in 
a tangent space approximation, and so do not have a finite measure space. 6 That is, there 
is a classical entropy S c (t) associated with some coarse-graining, that in the limit of small 
cell size and late times grows as 

S c (t) « h KS t = I diKi J t =>■ S c (nT) « 2 np , (69) 

\«l> 0 / 

which is exactly what we found in the quantum calculations, see Eqs. (33) and (48). 

We thus see that the asymptotic growth rate of the entanglement entropy is equal to 
the sum of the positive Lyapunov exponents, implying the asymptotic convergence of the 
classical entropy and entanglement entropy. 

VII. DISCUSSION 

In this paper we studied entanglement entropy growth under a Bogolubov transforma¬ 
tion that entangled two degrees of freedom and their canonical conjugates. To understand 
entanglement entropy we traced over half of the degrees of freedom. The rate growth of the 
entanglement entropy ended up being asymptotically equal to the sum of the two positive 
Lyapunov exponents of the system. This was true for all factorizations except on sets of 
measure zero. 

5 For more general systems, the relationship is given by the Ruelle inequality [42-44]. 

6 There are various rigorous definitions of local entropy and entropy for non-compact or non-finite measure 
spaces, e.g., [45, 46], but here we do not attempt to use these, and so essentially take Eq. (69) as our 
definition of S c . 
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We have seen that in the simple system we studied here there is a strong relationship 
between a notion of classical dynamical entropy and the growth of entanglement entropy. 
Namely, the former bounds the latter and they are almost always asymptotically equal. This 
is an indication that the classical dynamics, as encoded by dynamical systems quantities like 
the entropy, contains significant information about the associated quantum dynamics and 
the rates of scrambling for such systems, in line with many previous results on quantum 
dynamical systems. We note here though, that we considered almost-classical initial condi¬ 
tions. It would be interesting to analyze the dynamics of more intrinsically quantum states, 
like the superposition of widely separated coherent states. 

Our results can be generalized to more degrees of freedom, like the general linearized 
dynamics near periodic orbits for higher dimensional chaotic Hamiltonian systems. In the 
simplest cases, the degrees of freedom can be treated as independent and we expect that the 
growth of entanglement entropy will again asymptotically converge to the classical entropy. 
For example, if we partition the set of coordinates into two sets, one of k coordinates and 
one of n — k > k coordinates (where by coordinates here we mean canonical pairs written 
as ladder operator pairs), we expect the entanglement entropy for the factorization should 
grow asymptotically like the sum over the 2k largest positive Lyapunov exponents, except 
on sets of measure zero for the possible factorizations. This comes from the additivity of 
entanglement entropy for independent systems. 

Note that in these systems the Lyapunov exponents are paired: for every positive Lya¬ 
punov exponent k, there is a corresponding negative one —k, due to the conservation of the 
symplectic form. The maximum entropy growth would then be achieved when k is as close 
as possible to n — k, and we would get that the entropy growth would be the sum of all 
positive Lyapunov exponents. This is the same way that classical entropy (as described by 
Pesin’s formula) is supposed to behave for general dynamical systems. However, the details 
of the general case of N degrees of freedom remain for future work. 

There are other definitions of entropy in quantum dynamical systems, e.g., [47, 48], and 
it would be interesting to understand the relationship between these and the entanglement 
entropy studied here, for the kinds of systems we are considering. Additionally, the dynamics 
studied here was confined to a linearized analysis about some periodic orbit. To extend these 
results beyond this, we will have to confront the effects of folding and the highly non-trivial 
topology of chaotic domains in phase space [31-33]. 
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Appendix A: General Bogolubov transformation preserving charge 


We consider a phase space with two oscillators, whose ladder operator algebra is given 
by the operators a\a and b\b. We will assume that aJ has charge +1 and U has charge 
— 1 and that a discrete linear dynamics (generated by a unitary operator) on this system 
preserves the U(l) charge. The commutation relations are 


[a, a^] = 1 = [b, Y] 

(Al) 

a\ Y] = 0 = [a\ b] etc. 

(A2) 


The most general transformation that is linear and unitary and preserves the algebra is 
given by a Bogoliubov transformation to oscillators a\a and b\b. The most general such 
linear transformation compatible with the U(l) symmetry is given by 



Taking adjoints we have that 



/ B * 
l C* D* 



(A3) 


(A4) 


and unitarity implies that the commutation relations are not changed between the a, b and 
a, b operators. We need to verify those of equation (Al) and the one on the left of equation 
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(A2). Using the commutation relations we find that 

[a, a t ] = AA*[a, a f ] + BB*[b\ b\ = \A\ 2 - \B\ 2 = 1 . (A5) 

We can thus parametrize the solutions of this equation as A = cosh (p) exp (i6 + i/3), B = 
sinh(p) exp (icf> + i/3) where we have added a redundancy of the description with the angle (3. 
We also get that 


[a\ 6 f ] = AC* [a\ a} + BD* [b, 6 f ] = BD* - AC* = 0 . 


(A6) 


Thus 


D A* exp(— id) coshp 

C = ~B* 


(AT) 


exp(—20) sinh p 

Finally, we also get that \D\- — \C\ 2 = 1, so D, C differ from A*, B* by a phase. We choose 
the phase so that it is easily expressible in terms of f3, so that we get 


A B \ I e ie cosh p sinh p 

= exp (i(3) 


C D 


e ^ sinh p e ie cosh p 


(AS) 


Appendix B: Computing entropy for Gaussian density matrices 

Consider the holomorphic quantization of the harmonic oscillator, also called the Bargmann 
representation, see, e.g., [38, 49] for reviews. An arbitrary state vector |/) in the Hilbert 
space of the harmonic oscillator is represented by a holomorphic function f(z) and the inner 
product is given by 

(/Is) = 2 ~~ J d 2 £exp (~zz)f*{z)g(z) . (Bl) 

Here we define the integral by taking z = x + iy, z = x — iy, and the measure d 2 z = dz Adz = 
2idxAdy. The limits of integration for x, y are from —oo to oo. For example, take the states 
| z n ) represented by z n . With this inner product we have that 

z m \z n ) = —— [ d 2 zexp(—zz)z m z n = m\5 mn 
2m J 

The raising and lowering operators of the harmonic oscillator act as multiplication by z 
and derivatives with respect to z, respectively. Orthonormal energy eigenstates | n) are 
represented by z n /\fn\. A general bounded linear operator O can be written as a formal 
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sum of ket-bra combinations and thus represented by an integral kernel O —>■ 0(z, w) whose 
action on |/) is defined by 



(B3) 


which produces a holomorphic function of z. The trace of an operator is given by 
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J d 2 z exp(—zz)0(z, z) . 


(B4) 


We assume that 0(z,w) is holomorphic in z and w, and hence can be written 


OO 


0(z,w) - y O nm z n w m . 


(B5) 


n,m =0 


The coefficients of this expansion are the matrix elements of O in the energy eigenbasis: 


O 


nm 


(n\0\m) . 


(B6) 


Determining the kernel of a given operator, and vice versa, is thus straightforward in princi¬ 


ple, but can be difficult in practice if the matrix elements are not simple functions of n and 


m. However, for Gaussian operators the problem is tractable because they can be specified 
by just a few parameters. 

We are specifically interested in density operators whose kernels have a Gaussian form: 


p(z, z) = B 0 exp (Biz 2 + B 2 zz + B{z 2 ) , 


(B7) 


which arose in our Eqs. (41) and (42) (here we treat z and z as independent). These density 
operators represent generalized squeezed states. To compute the entropy of these density 
operators we would like to find an exponential form so that we can write down log p. Their 
entropy was computed this way in [38] and we reproduce their computation here, giving 
more details of the derivation. To this end first consider the parametrization 



(B8) 


The parameters here correspond to the expectation values of certain simple operators: 




(BIO) 


(B9) 
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Normalization of p, Tr(p) = 1, implies that N = y 2 — |a;| 2 , while hermiticity requires y to be 
real and y > 0. Comparing with equation (B7) we find that 


— = B 1 

2 N 

(B12) 

i- — = b 2 

N 2 

(B13) 

(VN)- 1 = B 0 

(B14) 


We can easily solve for N, finding that 

N= [( j B 2 -1) 2 -4| j B 1 | 2 ]- 1 (B15) 


from which x, y readily follow. It will be useful below to define the following quantity 

2 _ f 1\ 2 | ,2 _ B 2 1 _ 41-BjJ 2 — (1 + B 2 ) 2 

X V 2) ^ B 2 + A 16\Bi\ 2 — 4(1 — B 2 ) 2 

The von Neumann entropy for a density operator is defined as 


(B16) 


s(p) = —TV(plogp). 


(B17) 


To compute logp we look for an explicitly exponential form for p (as an operator, not just 
its representation). The Gaussian form of the density operator representation in Eq. (B8) 
suggests a Gaussian ansatz for the operator itself, i.e., 


P 


C exp 



(^A(ciq) T o)cl) T Bo, T Bo) ) 


(BIS) 


We take A e M and A > \B\ to ensure convergence. Our starting assumption is that 
an operator of this form does indeed correspond to Eq. (B8). If so, then it is completely 
characterized by the expectation values of the operators cda, cd 2 and a 2 , so we will compute 
these and then perform the matching via Eqs. (B9) - (Bll). 

To compute these expectation values for the density operator in Eq. (B18), first notice 
that the form of the exponent can be simplified, since it is of the form of a standard thermal 
density operator after a Bogoliubov transfomation. We can use the inverse transformation 
to bring the exponent in Eq. (B18) to the standard form. Specifically, we write an inverse 
Bogoliubov transformation as 



cosh 9 

—e~^ sinh 9 


sinh 9 
cosh 9 



(B19) 
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Then, by taking d and (j) such that 


A = 2D cosh 29 
B = 2De~ i( ^ sinh 29 , 


(B20) 

(B21) 


where D = \\J A? — \B\ 2 , we have that 


p = C exp 


— - (A cosh 29 — \B \ sinh 2d) ( bb + + b^b) 


= C exp [-£>(66+ + b ] b)\ (B22) 


In this form, we can calculate expectation values of various operators. We start with 

(a 2 ) = Tr(a 2 p) = /cosh 2 9 b 2 — -e 1 ^ sinh 2d(6 +b + btf) + e 2u ^ sinh 2 9 D 2 \ . (B23) 


First, we notice that 


(b 2 ) = (& t2 > = 0, 


(B24) 


which one can see by inserting b 2 or lC 2 into the trace, expanding in a basis of eigenvectors 
of the number operator N b = tfb, and using orthogonality. Next, from Eq. (B22) we have 


p = Ce -D(btf+b'b) = Ce -D(2bib+1) = Ce -D e -2DN b _ 


(B25) 


First, we have that (1) = Trp = (C/2)(sinh D) 1 , which means that to normalize we should 
take C = 2 sinh D. Next, we can compute 


(Db) = Tr (tfbp) = Tr (bptf) = Ce D ^^(n\be 2DNb C\n) 

n =0 
oo 

= Ce~ D ^(n + 1 )(n+ l\e~ 2DNb \n + 1) 


77—0 

OO 


Ce D ^(n + 1)< 


—2D(n+\) 


n=0 


= Ce~ D e~ 2D 


,-2 D 


(1 — e~ 2D ) 2 


+ e 


-2D 


1 — e~ 2D 


C e 


-D 


4 sinh 2 D 
1 e~ D 


2 sinh D 


So we have 




(fttfo + bb t) = 2(6+6) + (1) = —— + 1 = cothD. 

sinh D 


(B26) 

(B27) 

(B28) 

(B29) 

(B30) 

(B31) 

(B32) 










From here we can match to Eqs. (B9) - (Bll): 


1 B 

(a 1 2 ) = — sinh 20 coth D = — — coth D = x (B33) 

and 

(aa)) = cosh 2 9(bb 1) + sinh 2 6{tfb) = cothD + ^ = y . (B34) 

Returning to the computation of the entropy, we have that 

S(p) = -(log p) = -(log C) + ^(aa f + a)a) + ^(a 2 ) + ^(a t2 ) • (B35) 

We ultimately want S(p) in terms of the original parameters Rj, which we will do by first 
writing things in terms of 

X = ^ coth D = \j{y~ t ;) 2 - M 2 , (B36) 

from which one can then use Eq. (B16). First, using hyperbolic trigonometric identities, 
one can show that — log C = — log(2 sinh D) = \ log(y + 4) -f 1 log(% — \). Then, solve the 
final equations in Eqs. (B33) and (B34) for B and A, respectively, and substitute them into 
Eq. (B35), along with the expressions for the expectation values in terms of x and y. One 


ends up with 

S{p) = i log(x + f) + \ log(x - I) + (^(y - ^) 2 - |x| 2 ^ (B37) 

= \ log(x + |) + \ log(x - |) + D coth D (B38) 

= \ log(x +\) + \ log(x - \) + X log () (B39) 

\X 2 / 

= (x + I) l°g(x + \)-(x-\) log(x - |) , (B40) 


where the penultimate equality also comes from hyperbolic trigonometric identities. This 
gives Eq. (46). 
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